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Abstract 

We describe an approach based on topology optimization that enables automatic discovery of wavelength- 
scale photonic structures for achieving high-efficiency second-harmonic generation (SHG). A key distinc¬ 
tion from previous formulation and designs that seek to maximize Purcell factors at individual frequencies 
is that our method not only aims to achieve frequency matching (across an entire octave) and large radiative 
lifetimes, but also optimizes the equally important nonlinear-coupling figure of merit j3, involving a compli¬ 
cated spatial overlap-integral between modes. We apply this method to the particular problem of optimizing 
micropost and grating-slab cavities (one-dimensional multilayered structures) and demonstrate that a vari¬ 
ety of material platforms can support modes with the requisite frequencies, large lifetimes Q > 10 4 , small 
modal volumes ~ (A/n) 3 , and extremely large B > HP 2 , leading to orders of magnitude enhancements 
in SHG efficiency compared to state of the art photonic designs. Such giant j3 alleviate the need for ultra¬ 
narrow linewidths and thus pave the way for wavelength-scale SHG devices with faster operating timescales 
and higher tolerance to fabrication imperfections. 

PACS numbers: Valid PACS appear here 
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Introduction .—Nonlinear optical processes mediated by second-order (x (2> ) nonlinearities play 
a crucial role in many photonic applications, including ultra-short pulse shaping [it 2], spec¬ 
troscopy [3], generation of novel frequencies and states of light [4-fj], and quantum informa¬ 
tion processing [7-[9]. Because nonlinearities are generally weak in bulk media, a well-known 
approach for lowering the power requirements of devices is to enhance nonlinear interactions 
by employing optical resonators that confine light for long times (dimensionless lifetimes Q ) in 
small volumes V llOKlMl. Microcavity resonators designed for on-chip, infrared applications offer 
some of the smallest confinement factors available, but their implementation in practical devices 
has been largely hampered by the difficult task of identifying wavelength-scale (V ~ A 3 ) struc¬ 
tures supporting long-lived, resonant modes at widely separated wavelengths and satisfying rigid 


frequency-matching and mode-overlap constraints 115 . 


200 . 


In this letter, we extend a recently proposed formulation for the scalable topology optimization 
of microcavities, where every pixel of the geometry is a degree of freedom, to the problem of 
designing wavelength-scale photonic structures for second harmonic generation (SHG). We apply 
this approach to obtain novel micropost and grating microcavity designs supporting strongly cou¬ 
pled fundamental and harmonic modes at infrared and visible wavelengths with relatively large 
lifetimes Qi,Q 2 > 10 4 . In contrast to recently proposed designs based on known, linear cavity 
structures hand-tailored to maximize the Purcell factors or mode volumes of individual resonances, 


e.g. ring resonators [J_7 


21 


23] and nanobeam cavities [19, 


24], our designs ensure frequency 


matching and small confinement factors while also simultaneously maximizing the SHG enhance¬ 
ment factor Q 2 1 Q 2 \P\ 2 to yield orders of magnitude improvements in the nonlinear coupling 3 
described by ([3]) and determined by a special overlap integral between the modes. These particular 
optimizations of multilayer stacks illustrate the benefits of our formalism in an approachable and 
experimentally feasible setting, laying the framework for future topology optimization of 2D/3D 
slab structures that are sure to yield even further improvements. 


Structure 

h x x h y x h z (Af) 

A (/rm) 

(Qi, Qd 

(Qi ad , Qt) 

P 

FOMi 

fom 2 

(1) AlGaAs/ALCL micropost 

8.4 x 3.5 x 0.84 

1.5-0.75 

(5000, 1000) 

(1.4 x 10 5 ,1.3 x 10 5 ) 

0.018 

7.5 x 10 6 

8.3 x 10 11 

(2) GaAs gratings in Si0 2 

5.4 x 3.5 x 0.60 

OO 

1 

o 

so 

(5000, 1000) 

(5.2 x 10 4 , 7100) 

0.020 

7 x 10 6 

7.5 x 10® 

(3) LN gratings in air 

5.4 x 3.5 x 0.80 

o 

bo 

1 

o 

4 ^ 

(5000, 1000) 

(6700, 2400) 

0.030 

8.4 x 10 5 

9.7 x 10 7 


TABLE I. SHG figures of merit for topology-optimized micropost and grating cavities of different material 
systems. 
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Most experimental demonstrations of SHG in chip-based photonic systems Ilia. 
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operate in the so-called small-signal regime of weak nonlinearities, where the lack of pump deple¬ 
tion leads to the well-known quadratic scaling of harmonic output power with incident power [33]. 
In situations involving all-resonant conversion, where confinement and long interaction times lead 


to strong nonlinearities and non-negligible down conversion 11121. 


20|], the maximum achievable 


conversion efficiency r] = -y^r 
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occurs at a critical input power [20Q, 
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where Xeir' s the effective nonlinear susceptibility of the medium [SM], 0 = 
dimensionless quality factor (ignoring material absorption) incorporating radiative decay and 
coupling to an input/output channel The dimensionless coupling coefficient 3 is given by a 
complicated, spatial overlap-integral involving the fundamental and harmonic modes [SM], 

J dr e{r)E; El 
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where e(r) = 1 inside the nonlinear medium and zero elsewhere. Based on the above expressions 
one can define the following dimensionless figures of merit, 
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where FOMi represents the e. 


ficiency per power, often quoted in the so-called undepleted regime 


of low-power conversion [33y, and FOM 2 represents limits to power enhancement. Note that 


for a given radiative loss rate, FOMi is maximized when the modes are critically coupled, Q = 


Q n 


2 , with the absolute maximum occurring in the absence of radiative losses, g rad —> oo, or 
equivalently, when FOM 2 is maximized. From either FOM, it is clear that apart from frequency 
matching and lifetime engineering, the design of optimal SHG cavities rests on achieving a large 
nonlinear coupling /3. 

Optimal designs .— Table I characterizes the FOMs of some of our newly discovered microcav¬ 
ity designs, involving simple micropost and gratings structures of various x^ materials, including 
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GaAs, AlGaAs and LiNb 03 . The low-index material layers of the microposts consist of alumina 
(AI 2 O 3 ), while gratings are embedded in either silica or air (see supplement for detailed specifi¬ 
cations). Note that in addition to their performance characteristics, these structures are also sig¬ 
nificantly different from those obtained by conventional methods in that traditional designs often 
involve rings tfl7L Il8i| . periodic structures or tapered defects [24], which tend to ignore or sacri¬ 
fice /3 in favor of increased lifetimes and for which it is also difficult to obtain widely separated 
modes 1119(1 . Figure Q] illustrates one of the optimized structures—a doubly-resonant rectangular 
micropost cavity with alternating AlGaAs/AI 2 O 3 layers—along with spatial profiles of the funda¬ 
mental and harmonic modes. It differs from conventional microposts in that it does not consist of 
periodic bi-layers yet it supports two localized modes at precisely Ai = 1.5 /rm and A 2 = Ai/2. In 
addition to having large Q rad > 10 5 and small V ~ (\\/n) ', the structure exhibits an ultra-large 
nonlinear coupling j3 ~ 0.018 that is almost an order of magnitude larger than the best overlap 
found in the literature (see Fig. [2]). From an experimental point of view, the micropost system is of 
particular interest because it can be realized by a combination of existing fabrication techniques 
such as molecular beam epitaxy, atomic layer deposition, selective oxidation and electron-beam 
lithography [25]. Additionally, the micropost cavity can be naturally integrated with quantum dots 
and quantum wells for cavity QED applications [260. Similar to other wavelength-scale struc¬ 
tures, the operational bandwidths of these structures are limited by radiative losses in the lateral 


direction lllOl. 


25, 


27], but their ultra-large overlap factors more than compensate for the increased 


bandwidth, which ultimately may prove beneficial in experiments subject to fabrication imperfec¬ 
tions and for large-bandwidth applications [ 1 , 2 , 6 , 28]. 


To understand the mechanism of improvement in /3, it is instructive to consider the spatial pro¬ 
files of interacting modes. Figure [TJ? plots the //-components of the electric fields in the .x;j-planc 
against the background structure. Since j3 is a net total of positive and negative contributions 
coming from the local overlap factor EfE 2 in the presence of nonlinearity, not all local contribu¬ 
tions are useful for SHG conversion. Most notably, one observes that the positions of negative 
anti-nodes of E 2 (light red regions) coincide with either the nodes of E } or alumina layers (where 
X (2) = 0), minimizing negative contributions to the integrated overlap. In other words, improve¬ 
ments in j3 do not arise purely due to tight modal confinement but also from the constructive 
overlap of the modes enabled by the strategic positioning of field extrema along the structure. 

Based on the tabulated FOMs (Table I), the efficiencies and power requirements of realistic 
devices can be directly calculated. For example, assuming Xeff (AlGaAs) ~ 100 pnr/V [18], 
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FIG. 1. (a) Schematic illustration of a topology-optimized micropost cavity with alternating AlGaAs/Al 2 03 
layers and dimensions h x x h y x h z = 8.4 x 3.5 x 0.84 (Af). For detailed structural specifications, please 
refer to the supplement, (b) x-z cross-section of the E y components of two localized modes of frequencies 
Ai = 1.5pm and A 9 = Ai/2 and large spatial overlap E* 2 El 


the AlGaAs/Al 2 03 micropost cavity (Fig. Q} yields an efficiency of ^0r~ = 2.7 x 10 4 /W in 
the undepleted regime when the modes are critically coupled, Q = For larger operational 
bandwidths, e.g. Qi = 5000 and Q 2 = 1000, we find that P pj> ut = 16/W. When the system is 
in the depleted regime and critically coupled, we find that a maximum efficiency of 25% can be 
achieved at P{ ui m 0.15 rnW whereas assuming smaller Q 1 = 5000 and () 2 = 1000, a maximum 
efficiency of 96% can be achieved at P" 1 * 0.96 W. 


Comparison against previous designs .— Table II summarizes various performance character¬ 
istics, including the aforementioned FOM, for a handful of previously studied geometries with 
length-scales spanning from mm to a few wavelengths (microns). Fig 2 demonstrates a trend 
among these geometries towards increasing ft and decreasing Q rad as device sizes decrease. Max¬ 
imizing /? in millimeter-to-centimeter scale bulky media translates to the well-known problem of 


phase-matching the momenta or propagation constants of the modes [33]. In this category, tra¬ 
ditional WGMRs offer a viable platform for achieving high-efficiency conversion L29D; however, 
their ultra-large lifetimes (critically dependent upon material-specific polishing techniques), large 
sizes (millimeter length-scales), and extremely weak nonlinear coupling (large mode volumes) 
render them far-from optimal chip-scale devices. Although miniature WGMRs such as microdisk 


and microring resonators 117 . 


30, 


320 show increased promise due to their smaller mode vol- 
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FIG. 2. Scatter plot of (Q| ad )" Q™ d versus nonlinear overlap |/3| 2 for representative geometries, including 


WGMRs [29], micro- and nano-ring resonators [ 17 . 


18 


, photonic crystal slab and nanobeam cavities 11Q . 


19TI . A trend towards decreasing lifetimes and increasing overlaps is readily observed as devices become 
increasingly smaller. Meanwhile, it remains an open problem to discover structures with high Q s, small V s 
and large |/3| (shaded region). 


umes, improvements in (3 are still hardly sufficient for achieving high efficiencies at low powers. 
Ultra-compact nanophotonic resonators such as the recently proposed nanorings [180, 2D pho- 


.... 


o 


tonic crystal defects 111611 , and nanobeam cavities [19], possess even smaller mode volumes but 
prove challenging for design due to the difficulty of finding well-confined modes at both the fun¬ 
damental and second harmonic frequencies [160. Even when two such resonances can be found 
by fine-tuning a limited set of geometric parameters [ 18l 119], the frequency-matching constraint 
invariably leads to sub-optimal spatial overlaps which severely limits the maximal achievable /?. 

Comparing Tables I and II, one observes that for a comparable Q, the topology-optimized 
structures perform significantly better in both FOMj and FOM 2 than any conventional geometry, 
with the exception of the LN gratings, whose low () rml lead to slightly lower FOM 2 . Generally, the 
optimized microposts and gratings perform better by virtue of a large and robust ft which, notably, 
is significantly larger than that of existing designs. Here, we have not included in our comparison 


those structures which ac 


lieve non-negligible SHG by special poling techniques and/or quasi¬ 


phase matching methods [32-34], though their performance is still sub-optimal compared to the 


topology-optimized designs. Such methods are highly material-dependent and are thus not readily 
applicable to other material platforms; instead, ours is a purely geometrical topology optimization 
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FIG. 3. Work flow of the design process. The degrees of freedom in our problem consist of all the pixels 
along x-direction in a 2D computational domain. Starting from the vacuum or a uniform slab, the optimiza¬ 
tion seeks to develop an optimal pattern of material layers (with a fixed thickness in the z-direction) that can 
tightly confine light at the desired frequencies while ensuring maximal spatial overlap between the confined 
modes. The developed 2D cross-sectional patterns is truncated at a finite width in the y-dircction to produce 
a fully three-dimensional micropost or grating cavity which is then simulated by FDTD methods to extract 
the resonant frequencies, quality factors, eigenmodes and corresponding modal overlaps. Here, it must be 
emphasized that we merely performed one-dimensional optimization (within a 2D computational problem) 
because of limited computational resources; consequently, our design space is severely constrained. 


technique applicable to any material system. 

Optimization formulation .— Optimization techniques have been regularly employed by the 
photonic device community, primarily for fine-tuning the characteristics of a pre-determined ge¬ 
ometry; the majority of these techniques involve probabilistic Monte-Carlo algorithms such as par¬ 


ticle swarms, simulated annealing and genetic algorithms [J35 


-37] - While some of these gradient- 


free methods have been used to uncover a few unexpected results out of a limited number of 
degrees of freedom (DOF) [38], gradient-based topology optimization methods efficiently handle 
a far larger design space, typically considering every pixel or voxel as a DOF in an extensive 2D 
or 3D computational domain, giving rise to novel topologies and geometries that might have been 


difficult to conceive from conventional intuition alone. 


mization were primarily focused on mechanical problems[39Q and only recently have they 


he early applications of topology opti- 

)ecn ex¬ 


panded to consider photonic systems, though largely limited to linear device designs [38 


4Q-44]. 
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Structure 

A (/im) 

(Qi, Q 2 ) 

(Qf d , Qt) 

P 

FOMi 

fom 2 

LN WGM resonator f291 

1.064 - 0.532 

(3.4 x 10 7 , -) 

(6.8 x 10 7 ,-) 

- 

~10 10 

- 

AIN microring U71 

1.55 - 0.775 

(~ 10 4 , — 5000) 

- 

- 

2.6 x 10 5 

- 

GaP PhC slab [lfi] * 

1.485 - 0.742 

(« 6000, -) 

- 

- 

« 2 x 10 5 

- 

GaAs PhC nanobeam IT91 

1.7 — 0.91'*' 

(5000, 1000) 

(> 10 6 , 4000) 

0.00021 

820 

1.8 x 10 8 

1.8-0.91 

(5000, 1000) 

(6 x 10 4 , 4000) 

0.00012 

227 

2.1 x 10 5 

AlGaAs nanoring H81 

1.55 - 0.775 

(5000, 1000) 

(10 4 , > 10 6 ) 

0.004 

10 5 

1.6 x 10 9 


TABLE II. SHG figures of merit, including the frequencies A, overall and radiative quality factors Q. Q md 
and nonlinear coupling j3 of the fundamental and harmonic modes, of representative geometries. Also 
shown are the FOMi and FOM 2 figures of merit described in ©, ([5]). 


* SHG occurs between a localized defect mode (at the fundamental frequency) andh an extended index- 
guided mode of the PhC. 
i Resonant frequencies are mismatched. 


In what follows, we describe a technique for gradient-based topology optimization of nonlinear 
wavelength-scale frequency converters. 


Recent work [41] considered topology optimization of the cavity Purcell factor by exploiting 


the concept of local density of states (LDOS). In particular, this previous formulation exploited 
the equivalency between LDOS and the power radiated by a point dipole in order to reduce 
Purcell-factor maximization problems to a series of small scattering calculations. Defining the 
objective max £ - / (e(r);ce) = — Re[J rfrJ* • E], it follows that E can be found by solving the 
frequency domain Maxwell’s equations A4E = icoJ, where A4 is the Maxwell operator [SM] and 
J = 5 (r — r 0 ) etj. The maximization is then performed over a finely discretized space defined 
by the normalized dielectric function, {e a = e(r a ), a o- (iAx,jAy,kAz)}. A key realiza¬ 
tion in [4 lj] is that instead of maximizing the LDOS at a single discrete frequency u, a better- 
posed problem is that of maximizing the frequency-averaged / in the vicinity of to, denoted by 
(/) = J doj' W(u/; u;, T)f(u'), where W is a weight function defined over some specified band¬ 
width T. Using contour integration techniques, the frequency integral can be conveniently replaced 
by a single evaluation of / at a complex frequency to + iT [41]. For a fixed T, the frequency aver¬ 
age effectively shifts the algorithm in favor of minimizing V over maximizing 0\ the latter can be 
enhanced over the course of the optimization by gradually winding down the averaging bandwidth 
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r [41|]. A major merit of the frequency-averaged LDOS formulation is that it features a mathe¬ 
matically well-posed objective as opposed to a direct maximization of the cavity Purcell factor y, 
allowing for rapid convergence of the optimization algorithm into an extremal solution. 


As suggested in [41], a simple extension of the optimization problem from single- to multi- 


mode cavities maximizes the minimum of a collection of LDOS at different frequencies. Applying 
such an approach to the problem of SHG, the optimization objective becomes: max^min LDOS(u;i), LDOS(2cci' 
which would require solving two separate scattering problems, M ]Ei = Ji and A4 2 E 2 = J 2 , 
for the two distinct point sources Ji, J 2 at ui\ and cc 2 = 2cci respectively. However, as discussed 
before, rather than maximizing the Purcell factor at individual resonances, the key to realizing 
optimal SHG is to maximize the overlap integral /? between Ei and E 2 , described by ([3]). Here, 
we suggest an elegant way to incorporate j3 by coupling the two scattering problems. In particu¬ 
lar, we consider not a point dipole but an extended source J 2 ~ Ef at cc 2 and optimize a single 


instead of two otherwise unrelated LDOS. 


combined radiated power / = — Re f dr ■ E 2 
The advantage of this approach is that / yields precisely the (3 parameter along with any resonant 
enhancement factors (~ Q/V) in Ei and E 2 . Intuitively, J 2 can be thought of as a nonlinear 
polarization current induced by Ei in the presence of the second order susceptibility tensor X (2 '3 


and in particular is given by J 2i = e(r) Yljk xifkE Xj E xk where the indices i, j, k run over the 


Cartesian coordinates. In general, m i xes polarizations and hence / is a sum of different con¬ 
tributions from various polarization-combinations. In what follows and for simplicity, we focus 
on the simplest case in which Ei and E 2 have the same polarization, corresponding to a diagonal 
X i2) tensor determined by a scalar x^a ■ Such an arrangement can be obtained for example by 
proper alignment of the crystal orientation axes [SM]. With this simplification, the generalization 
of the linear topology-optimization problem to the case of SHG becomes: 


max ?Q (/(e a ;o;i)) = -Re 


J* 2 • E 2 dr^ , 


AtiEi — iwiJi, 

A4 2 E 2 = icc 2 J 2 , cc 2 = 2cci 


( 6 ) 


9 









where 


Ji = 5(r Q - r 0 )ej, j e {x,y,z} 

J2 = e(r a )^-ej, 

Mi = Vx-Vx - ei(r a )oji, l = 1,2 

/i 

^ 111 T (^d l £ m ) ■ G [0, 1], 


and where e d denotes the dielectric contrast of the nonlinear medium and e m is that of the sur¬ 
rounding linear medium. Note that e a is allowed to vary continuously between 0 and 1 whereas 
the intermediate values can be penalized by so-called threshold projection filters [45fl. The scatter¬ 
ing framework makes it straightforward to calculate the derivatives off (and possible functional 
constraints) with respective to e a via the adjoint variable method [39-4lfl. The optimization prob¬ 
lem can then be solved by any of the many powerful algorithms for convcx, con servat i vc, separable 
approximations, such as the well-known method of moving asymptotes [460. 

For computational convenience, the optimization is carried out using a 2D computational cell 
(in the xz- plane), though the resulting optimized structures are given a finite transverse extension 
h y (along the y direction) to make realistic 3D devices (see Fig. [3]). In principle, the wider the 
transverse dimension, the better the cavity quality factors since they are closer to their 2D limit 
which only consists of radiation loss in the z direction; however, as h y increases, B decreases 
due to increasing mode volumes. In practice, we chose h y on the order of a few vacuum wave¬ 
lengths so as not to greatly compromise either Q or B. We then analyze the 3D structures via 
rigorous FDTD simulations to determine the resonant lifetimes and modal overlaps. By virtue 
of our optimization scheme, we invariably find that frequency matching is satisfied to within the 
mode linewidths. We note that our optimization method seeks to maximize the intrinsic geometric 
parameters such as Q md and B of an un-loaded cavity whereas the loaded cavity lifetime Q de¬ 
pends on the choice of coupling mechanism, e.g. free-space, fiber, or waveguide coupling, and is 
therefore an external parameter that can be considered independently of the optimization. When 
evaluating the performance characteristics such as FOM|, we assume total operational lifetimes 
Q 1 = 5000, Q 2 = 1000. In the optimized structures, it is interesting to note the appearance of 
deeply sub-wavelength features ~ 1 — 5% of —, creating a kind of metamaterial in the optimiza¬ 
tion direction; these arise during the optimization process regardless of starting conditions due to 
the low-dimensionality of the problem. We find that these features are not easily removable as 
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their absence greatly perturbs the quality factors and frequency matching. 


Concluding remarks .— We have presented a formulation that allows for large-scale optimiza¬ 
tion of SHG. Applied to simple micropost and grating structures, our approach yields new classes 
of microcavities with stronger performance metrics over existing designs. One potentially chal¬ 
lenging aspect for fabrication in the case of gratings is the presence of deeply sub-wavelength 
features, which would require difficult high-aspect-ratio etching or growth techniques. This is not 
an issue for the micropost cavities since each material layer can be grown/deposited to a nearly ar¬ 
bitrary thickness II25L 12611. Another caveat about wavelength-scale cavities is that they are sensitive 
to structural perturbations near the cavity center, where most of the field resides. In our optimized 
structures, the figures of merit are robust to within ~ ±20 nm variations (approximately one com¬ 
putational pixel). One possible way to constrain the optimization to ensure some minimum spatial 


feature and ro 


justness is to exploit so-called regularization filters and worst-case optimization 


techniques [45], which we will consider in future work. 


Our results provide just a glimpse of the kinds of designs that could be realized in structures 
with 2D and 3D variations, where we expect even better performance metrics due to the larger de¬ 
sign space. In fact, preliminary application of our formulation to 2D systems reveals overlap fac¬ 
tors and lifetimes at least one order of magnitude larger than those attained here. Apart from SHG 
optimization, our approach can be generalized to consider other nonlinear processes, even those 
involving more than two frequencies [SM]. Preliminary investigations reveal orders-of-magnitude 
improvements in the efficiency of third harmonic and sum-frequency generation processes. These 
findings, together with higher-dimensional applications, will be presented in future work. 


* zlin@seas.harvard.edu 
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Abstract 

We review the temporal-coupled mode equations describing second harmonic generation in doubly res¬ 
onant cavities and motivate the dimensionless nonlinear coupling j3 described in Eq. 3 of the main text. 
We provide further details on the topology optimization formulation for second harmonic generation and 
describe generalizations to other nonlinear processes. Finally, we present more detailed descriptions of the 
optimized nricropost and gratings cavities. 
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COUPLED-MODE THEORY FOR SECOND HARMONIC GENERATION 


The temporal coupled mode equations describing second harmonic generation (SHG) in a 
doubly-resonant cavity coupled to a channel are [1]: 
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such that |dfc| is the modal energy in the cavity and \Sk±\ is the input/output power in the waveg¬ 
uide, and where Qk and Q” d denote the total and radiative quality factors corresponding to mode 
k. The nonlinear coupling coefficient Bi, obtained from perturbation theory [11], is given by: 
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with B 2 = d\l‘l far off from material resonances where Kleinman symmetry is valid [2], as re¬ 
quired by conservation of energy [1]. In general, the overlap integral in the numerator is a sum 
of products between different E'-field polarizations weighted by off-diagonal components of the 
nonlinear y (2) tensor. For simplicity, however, in the main text we only consider the simple case 
of diagonal x m involving same-polarization interactions described by an effective x^ 2 J , resulting 
from an appropriate orientation of the crystal axes of the nonlinear material. All of these consid¬ 
erations suggest a simple dimensionless normalization of B, given by: 
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such that /? 2 = 4 Bxfn /\AoAf. As defined in the text, e(r) = 1 for nonlinear dielectric and 
e(r) = 0 for the surrounding linear medium. 

Most SHG experiments operate in the small-signal regime of small input powers, leading to 
negligible down-conversion and pump depletion. Ignoring the down-conversion or /3i term in 
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Eq.[Q one obtain the following simple expression for the second harmonic output power: 
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In the limit of large up-conversion and non-negligible down-conversion, solution of the coupled¬ 
mode equations yields the maximum efficiency (defined as 77 = P 2 ai / Pj n ) and corresponding 
critical power [lfl: 
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FORMULATION FOR TOPOLOGY OPTIMIZATION OF ARBITRARY NONLINEAR FRE¬ 
QUENCY CONVERSION PROCESS 


Nonlinear frequency conversion processes can be viewed as frequency mixing schemes in 
which two or more constituent photons at a set of frequencies {ui n } interact to produce an out¬ 
put photon at frequency Q such that Q = Y^ n c n uj n , where the photon number coefficients {c n } 
can be either negative or positive, depending on whether the corresponding photons are created 
or destroyed in the process. In general, because the optical nonlinear response of materials is a 
tensor and hence the frequency conversion process mixes different polarizations [2]. However, for 
notational simplicity, we will describe our optimization problem only for a single component of 
the susceptibility tensor. If one wishes to consider the full tensorial properties, one can easily add 
extra optimization terms (weighted by the tensor components) by following the same approach 
described below. 

Given a specific nonlinear tensor component Xijk..., where i,j, k, ... E {x, y, z}, mediating an 
interaction between the polarization components E^Q) and E^, E 2 k, •••, we begin with a collec¬ 
tion of point dipole currents, each at the constituent frequency ui n , n E { 1 , 2 ,...} and positioned 
at the center of the computational cell r', such that J n = e nu 5(r — r'), where e nv G {e^, e 2 fc,...} 
is a polarization vector chosen so as to excite the desired f?-field polarization components of the 
corresponding mode. Given the choice of incident currents J n , we solve Maxwell’s equations to 
obtain the corresponding constituent electric-field response E n , from which one can construct a 


nonlinear polarization current J(f2) = e(r) f] n Eni^ ( * ) e l , where E nv = E 


• e r 


and J( 0 ) has 
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polarization e, generally different from the constituent polarizations e nu . Here, (*) denotes com¬ 
plex conjugation for negative c n and no conjugation otherwise. Finally, we maximize the radiated 


power —Re 


f J(fi)* • E(fi) dr from J(0). 


The formulation is now given by: 
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In practice, we maximize a frequency-averaged version of the power output (f(u)) rather than 
f(u) itself since the latter has poor convergence [30, i.e., we maximize (/) = J du' W(o/; u, r)/(u/), 
where we simply choose the weighting function W to be a simple lorentzian with the desired res¬ 
onance u) and a certain linewidth I\ 


W(a/) = 


t/tt 


( 10 ) 
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(Note that this linewidth is not necessarily the same as the intrinsic radiative linewidth of the 
cavity in that T is a computational artifice introduced to aid rapid convergence [30.) In Ref. [30, 
it is shown by means of contour integration that this averaging is equivalent to evaluating / at a 
complex frequency f(ui + iT), also equivalent to adding a uniform loss iT/ui to /i(r) and e(r). In 
our implementation of the optimization process, we typically begin with a large T which affords 
rapid convergence to a stable geometry in a few hundred iterations. T is then decreased by an order 
of magnitude every time the optimization converges until T ~ 10 -5 at which point the structure 
settles to within a linewidth T of the desired frequencies (perfect frequency matching). 

Application of this formulation to the problem of second harmonic generation is straight¬ 
forward and described in the main text, in which case Q = u 2 = 2ui and J(H) = J 2 = 
e(r) (ey ■ Ei) 2 e 2 j. Note that for the structures described in the text, we chose e 1? = e 2l = e y . In 
addition to the problem statement of Eq. [9] the optimization algorithm also benefits from gradient 
information of the objective function, which exploits the adjoint variable method [3—5]. Here, we 
simply quote the result for the gradient of our SHG objective function (dropping the polarization 
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where the functions Uk are solutions of the following scattering problems: 


AfiWi = eE*E 2 , 
M.2^2 — tEfi 
Aiiu 3 = 2eE^U2- 


OPTIMIZED ID CAVITY DESIGNS 


Figure Q] shows the dielectric and E y field profiles of the three optimized structures, includ¬ 
ing (a) AlGaAs/Al 2 0 3 micropost, (b) GaAs gratings in Si0 2 , and (3) LN gratings in air. Along 
the x cross-section, each computational pixel of thickness A represents either a high dielectric 
(nonlinear) or low dielectric (linear) material. For example, in the AlGaAs/Al 2 03 micropost cav¬ 
ity (assuming ri\ (AlGaAs) = 3.02 and 77, 2 (AlGaAs) = 3.18 for AlGaAs with 70% A1 [6], and 
n(Al 2 0 3 ) = 1.7), we took A = 0.015 
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(a) AlGaAs/Al 2 0 3 micropost 



x/A 1 




(b) GaAs gratings in glass 
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(c) LN gratings in air 




FIG. 1. Cross-sectional dielectric profiles and electric field distributions for AIGaAs/ALO.a micropost (a), 
GaAs gratings in Si02 (b) and LN gratings in air (c). 
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